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ABSTRACT 

We present 23 candidate X-ray binaries with luminosities down to 1.8 x 10'^^ erg s~^, in the nearby 
starburst galaxy NGC 4449, from observations totaling 105 ksec taken with the ACIS-S instrument on 
the Chandra Space Telescope. We determine count rates, luminosities, and colors for each source, and 
perform spectral fits for sources with sufficient counts. We also compile a new catalog of 129 compact 
star clusters in NGC 4449 from high resolution, multi-band optical images taken with the Hubble 
Space Telescope, doubling the number of clusters known in this galaxy. The UBVI,Ha luminosities 
of each cluster are compared with predictions from stellar evolution models to estimate their ages and 
masses. We find strong evidence for a population of very young massive, black-hole binaries, which 
comprise nearly 50% of the detected X-ray binaries in NGC 4449. Approximately a third of these 
remain within their parent star clusters, which formed r ^ 6 — 8 Myr ago, while others have likely been 
ejected from their parent clusters. We also find evidence for a population of somewhat older X-ray 
binaries, including both supergiant and Be-binaries, which appear to be associated with somewhat 
older T « 100—400 Myr star clusters, and one X-ray binary in an ancient (r w 10 Gyr) globular cluster. 
Our results suggest that detailed information on star clusters can significantly improve constraints on 
X-ray binary populations in star-forming galaxies. 

Subject headings: galaxies: individual (NGC 4449) — galaxies: star clusters — binaries: close — stars: 
evolution — X-rays 



1. INTRODUCTION 

Images of nearby starburst galaxies taken in X-rays 
with the Chandra Space Telescope are spectacular, show- 
ing a multitude of bright point sources. It is now gen- 
erally accepted that most of these sources are high mass 
X-ray binaries (HMXBs) produced during recent star for- 
mation. HMXBs are binaries where one member of the 
system is a compact object, either a black hole or neu- 
tron star, and the other is a young, massive star. X-ray 
emission is produced as material is accreted from the 
young "donor" star onto the compact object. HMXBs 
can be divided into two general categories based on the 
type of donor star: 1) a Be star (Be/X-ray binary), or 2) 
a supergiant star (SG/X-ray binary). 

Previous studies have suggested an intriguing connec- 
tion between HMXBs and young stellar clusters: many 
of the former are found close to, but not coincident with, 
young star clusters. These observations are broadly con- 
sistent with a scenario where X-ray binaries (XRBs) form 
in star clusters, but have sufficiently large ve locities that 
they a r e expelled from their parent cluster (jZezas et al.l 
(|20'02D . IKaaret et all (|2004) ). There are three mecha- 
nisms that might displace HMXBs from their birthsites. 
The binary could be given a "kick" during an asymmetric 
supernova explosion that forms the black hole or neutron 
star, or it could be ejected via dynamical inter actions 
with other stars in a dense cluster core (McSwa in et al.l 
12001 . A third scenario that could explain the observed 
spatial displacement between HMXBs and young star 
clusters is that X-ray binaries formed in clusters that 
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have since dispersed. 

In this paper, we investigate the relationship between 
HMXBs and star clusters using the best X-ray and op- 
tical observations currently available for the starburst 
galaxy NGC 4449. We analyse observations taken with 
the Chandra Space Telescope to select X-ray binaries in 
NGC 4449, and optical images taken with the Hubble 
Space Telescope {HST) to detect star clusters and to es- 
timate their ages and masses. An optical color image of 
NGC 4449 based on the HST images is shown in Fig- 
ure 1. Our primary goals are to constrain the nature 
of X-ray binaries and in particular HMXBs, and to ex- 
plore the relationship between HMXBs and star clusters 
in NGC 4449. 

The rest of this paper is organized as follows. In Sec- 
tion 2 we present the X-ray observations from Chan- 
dra, including data reduction, source detection and some 
basic models used for interpretation. In Section 3 we 
present the optical observations from the HST and a 
new catalog of compact star clusters, including their 
integrated colors and luminosities, and we derive their 
masses and ages. The mass and age distributions of the 
clusters and what they tell us about the formation and 
disruption of the clusters are the subject of Section 4, 
while Section 5 presents the spatial correlation between 
the clusters and the candidate HMXBs. Section 6 synthe- 
sizes these results to constrain the nature of the HMXBs 
in NGC 4449, and Section 7 summarizes the main results 
of this work. 

2. X-RAY OBSERVATIONS FROM CHANDRA 

NGC 4449 is an irregular, star-f orming galaxy fo cated 
at a distance of 3.82 ± 0.27 Mpc ()Annibali et al.| [2008'). 
With an integrated magnitude of Mb = —18.2 it is some- 
what more lumninous than the Large Magellanic Cloud 
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Figure 1. Color image of NGC 4449 produced using the H ST / ACS observations described in Section 3. The F435W (~ B band) image 
is shown in blue, the F555W (~ V band) image in green, and the red image combines the F814W (~ / band) and F658N (Hq) titers. 
North is up and east is to the left. This image is 5' along the long side. (Credit: NASA, ESA. A. Aloisi (STScI/ESA), and The Hubble 
Heritage (STScI/AURA)-ESA/Hubble Collaboration) 



(|Huiitej[l997l ). NGC 4449 has a current sta r formation 
rate of w 1.5 Mq yr~^ (jThronson et al.|[l987l ) and a near- 
solar present-da y gas abundance (12 + \og[0/H] — 8.83; 
iGrevesse fc Sau val (1998)). 

2.1. Data and Reduction 

We use three sets of archival Chandra observations 
of NGC 4449, with integration times of 30 ksec (Ob- 
sID: 2031, PI: Heckman), 15 ksec (ObsID: 10125. PI: 
Long), and 60 ksec (ObsID: 10875, PI: Long), to detect 
X-ray point sources in NGC 4449. Basic information for 
the three sets of observations is given in Table 1. The 
data were taken with the Advanced CCD Imaging Spec- 
trometer (ACIS) instrument on the Chandra telescope 
on February 4, 2001 in "faint" mode (ObsID: 2031), 
and March 4 and 7, 2010 in "vfaint" mode (ObsIDs: 
10125 and 10875). The galaxy was positioned on the 



back-illuminated S3 CCD chip. We processed the data 
using the Chandra Interactive Analysis of Observations 
(CIAO) software (version 4.2) and Chandra Calibration 
Data Base (CALDB) version 4.3.CQ, and restricted the 
data to the energy range between 0.3-8 keV. The obser- 
vations were filtered in three energy bands - 0.3-1 keV 
(soft), 1-2 keV (medium) and 2-8 keV (hard). 

2.2. Source Detection and X-ray Properties 

We use CIAO's Mexican-hat wavelet source detection 
routine wavdetect (iFreeman et al.|[2002 fl to create source 
lists. Wavelet scales of 1.4, 2, 4, 8, and 16 pixels and 
a detection threshold of 10~^ were used, which typically 
results in 1 spurious detection per million pixels. The 
output sources were examined visually to verify each de- 

^ http://cxc.harvard.edu/ciao/ 



XRBs and Star Clusters in NGC 4449 



3 



Table 1 

Chandra observations 



Obsid 


Date 


PI 


Exposure 
(ks) 


2031 


2001-02-04 


Heckman 


30 


10125 


2009-03-04 


Long 


15 


10875 


2009-03-07 


Long 


60 



tection, and to correct the source catalog when multiple 
detections occurred. A catalog of X-ray point sources 
detected in the Chandra observations is presented in Ta- 
ble [21 and the locations of these sources are shown in 
Figured 

A number of measurements and estimates are made 
for each source — total number of counts and the counts 
measured in the soft, medium, and hard energy bands. 
We calculate two X-ray colors, a "soft" color defined 
as HI = [M - S)/T, and a "hard" color defined as 
H2 = {H - M) /T, where 5, M and H are the total 
counts measured in th e soft, medium and hard bands 
(jPrestwich et al.l I2003D . and T is the total number of 
counts in all three bands. The luminosity of each source 
is estimated by fitting its spectrum with a power law 
model with a photon index T — 1.5. These values are 
compiled in Tabled 

There are eleven sources which have sufhcient counts 
(> 50) for a crude spectral fit. We were able to fit ten of 
these with simple X-ray spectral models, as summarized 
in Table [3] One source, X15, is in a region with high 
background from diffuse emission and we were unable to 
obtain a satisfactory fit. X-ray spectra and responses 
(including sensitivity of the instrument and CCD) were 
extracted using standard CIAO software. Spectra were 
grouped for a minimum of 15 counts per bin. The fits 
were performed using XSPEC V 12.0 over an energy 
range 0.2—8.0 keV. For each source, we tried three mod- 
els: a simple power law (PL), a MEKAL modeQ, and 
a multi-color disk (MCE0) model. A fixed foreground 
Galactic column nn = 1.5 x lO^'^ cm^^ was assumed in 
each fit. The fits allow for an additional variable column 
due to absorption intrinsic to the source. The best fit 
model and parameters for each source are given in Ta- 
ble d Observed X-ray fluxes and estimated luminosities 
were calculated for the 0.3 — 8 keV energy range (both 
uncorrected for absorption) using the best fit parameters 
in the table. 

2.3. X-ray Colors and Models 

After examining the resulting X-ray source catalog, 
we eliminate some sources from further consideration in 
this work for the following reasons. First, we remove 
all sources located beyond the optically luminous por- 
tion of NGC 4449 (R ~ 3-4 kpc), since these are al- 
most certainly background galaxies and not associated 
with NGC 4449 itself. We also eliminate one source 

^ A n emission spectrum from hot diffuse gas, e.g. ILiedahl et all 

^ A superposition of multi-temperature blackbody spectra ex- 
pected from optically-thick accretion disk; the m odel only con- 
strains the temperature of the inner disk Ti„, e.g. IMitsuda et al.l 



Table 2 

X-ray source catalogue 



ID 


RA 


DEC 


Soft 


Hard 










Color 


Color 


(erg s^"*") 


X8 


187.00308 


44.07569 


0.20 


0.00 


4.55E-I-36 


X14 


187.00572 


44.09145 


0.15 


-0.50 


9.41E-I-36 


X18 


187.01634 


44.09551 


0.01 


-0.09 


7.68E-I-37 


X13 


187.02855 


44.09120 


0.31 


-0.42 


9.21E-f36 


X6 


187.02988 


44.07074 


0.34 


0.14 


1.47E-I-37 


X29 


187.03043 


44.12264 


0.04 


-0.07 


3.34E-f37 


XIO 


187.03065 


44.08164 


-0.33 


-0.16 


5.67E-f37 


X12 


187.03891 


44.08566 


-0.06 


-0.32 


3.01E-I-38 


X20 


187.04055 


44.09808 


0.23 


-0.13 


1.01E-f38 


Xll 


187.04576 


44.08327 


0.34 


-0.26 


3.05E-I-36 


X23 


187.04671 


44.11066 


-0.24 


-0.26 


2.25E-I-37 


X15 


187.04985 


44.09188 


-0.44 


-0.04 


7.66E-f37 


X24 


187.04997 


44.11149 


0.02 


-0.29 


5.81E-f37 


X21 


187.05013 


44.09955 


-0.28 


-0.25 


1.19E-I-37 


X28 


187.05533 


44.11557 


-0.01 


-0.32 


6.93E-f37 


X31 


187.06829 


44.12765 


-0.08 


0.41 


9.53E-I-36 


X22 


187.07433 


44.10939 


0.41 


-0.13 


8.59E-I-38 


X19 


187.07917 


44.09585 


-0.90 


-0.02 


2.01E-I-37 


X41 


187.04336 


44.09946 


-1.00 


0.00 


1.43E-I-37 


X42 


187.06327 


44.08829 


-0.41 


-0.57 


1.78E-I-36 


X43 


187.08341 


44.10589 


0.41 


0.09 


3.36E-I-36 


X44 


187.05553 


44.11308 


0.13 


0.09 


1.13E-I-36 


X46 


187.04041 


44.08867 


-0.20 


-0.23 


7.59E-I-36 


X26^ 


187.04569 


44.11343 


0.11 


-0.35 


2.57E+38 



X-ray luminosities (0.3 — 8 kcV range) are derived by fitting a power 
law (r — 1.5) to the data. An assumed Galactic uh = 1-5 X 10^^ cm~^ 
is applied. 

^ Supernova remnant. 

that is a confirmed sup ernova remnant in NGC 4449 
(jPatnaude fc: Fesenll2003| ) , and another which appears to 
be a spurious detection of diffuse emission near the cen- 
ter of the galaxy. Two additional X-ray sources are co- 
incident with foreground stars and were also eliminated. 
Our final catalog contains 23 sources that we believe are 
X-ray binaries. These are listed in Tabled 

The observed HI vs. H2 colors of the X-ray binaries 
are shown in Figure [3] Two models are shown for refer- 
ence. The red line shows predictions from disk blackbody 
models with temperatures ranging from 0.1 — 1.0 keV. 
Black hole binaries typically have disk blackbody tem- 
peratures of kT ~ 1.0 keV. The orange line shows the 
effect of adding absorption to the disk blackbody mod- 
els. The green line represents a power law with increasing 
photon index from 1.0 — 3.0. Accreting low mass neutron 
star binaries typically have absorbed power law spectra 
with photon index ~ 2.0. Figure d suggests that NGC 
4449 has a mix of X-ray binary populations typical of a 
galaxy that has had on-going star formation, with some 
sources better following disk black-body models and oth- 
ers in the absorbed portion of the diagram. The nature 
of the X-ray binaries will be discussed in more detail in 
Section 6. 

3. OPTICAL OBSERVATIONS FROM HST 

It has been suggested that the dense environments 
found in compact star clusters may b e very efficient in 
producing XRBs fe.g.. iMcSwain et al.l (|2007D '). Regard- 
less of whether or not XRBs form within compact clus- 
ters, the clusters are good tracers of star formation in 
galaxies, and hence can provide important constraints 
on the XRB population. The goals of this section are 
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Figure 2. V band (F555W filter) image of NGC 4449 taken witli tlie ACS camera on the HST. North is up and east is to the left. The 
locations of star clusters selected from optical HST images are shown as circles. The locations of X-ray sources selected from Chandra 
X-ray observations are labeled and shown as crosses. 



two-fold: (1) select a new catalog of compact star clus- 
ters in NGC 4449, to correlate with the locations of the 
XRBs, and (2) measure luminosities and colors for any 
point source coincident with an XRB and hence likely 
to be the donor star (discussed in Sections 5.2 and 6.4). 
While a catalog of clusters in NG C 4449 already exists 
(jCelatt. Hunter fc Gallagheiil200lD . it is based on par- 
tial imaging taken with the WFPC2 instrument on board 
HST, and there are now deeper, higher resolution data 
with full coverage of NGC 4449 available from the Wide 
Field Channel (WFC) camera of the Advanced Camera 
for Surveys (ACS) instrument. 

3.1. Data and Photometry 

ACS/WFC imaging of two positions within NGC 4449 
was taken in the F435W (« B), F555W (« V band). 



F814W (« / band), and the F658N (Ha) fihers, on 
November 10-11, 2005 (Proposal ID GO: 10585, PI: 
Aloisi). Four individual exposures were taken in each 
filter at each pointing. The ACS/WFC has a pixel scale 
of 0."049, or a projected scale of 18.4 pc p er arcsecond at 
the a ssumed distance of 3.82 ± 0.27 Mpc (jAnnibali et al.l 
[2008f) to NGC 4449. We downloaded the ACS data from 
the Hubble Legacy Archive^ (HLA). The HLA combines 
the individual flatfielded exposures for a specific filter 
together using the PYRAF task Multidrizzle, and out- 
puts geometrically corrected images. For the NGC 4449 
observations, the HLA used stars from the U.S. Naval 
Observatory catalogue to astrometrically correct the im- 
ages. 

http://hla.stsci.edu/ 
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Table 3 

Model parameters for X-ray sources with > 50 counts 



ID 


Net 


Best Fit 






Photon 


T 1 


rp2 


/ d.o.f. 


Flux^ 






Counts 


Model 


( X 10 cm 


) 


index 


(IceV) 


(key) 




/ —1 — 2^ 
(erg s cm j 


(erg s 


X18 


141.3 


ABS*PL 


{4.5tt0)xl0- 


-2 


1.5±0.3 


- 


- 


2.4 / 7 


3.65E-14 


6.30E-I-37 


X29 


84.8 


ABS*PL 


(2.2+2^2^) X 10^ 


~2 


1 26+°'^ 






0.43 / 3 


2.48E-14 


4.30E-I-37 


XIO^' 


172.7 


ABS*(PL+MEKAL) 


{9.5lfo)xlO- 


-3 


2.1±0.4 




0.46l0;f 


7.0 / 7 


2.72E-14 


4.70E+37 


X12 


1197.7 


ABS*MCD 


(4.0lt5)xl0- 


-3 




0.60±0.5 




85.9 / 71 


1.58E-13 


2.76E+38 


X20 


201.4 


ABS*PL 


0.19±0.1 




1 7+0.4 






13.72 / 11 


4.80E-14 


8.35E+37 


X23*^ 


106.0 


ABS*MCD 








0.37±0.13 




5.28 / 4 


1.13E-14 


1.95E+37 


X24 


186.8 


ABS*PL 


10+006 










7.54 / 10 


3.20E-14 


5.66E+37 


X21 


58.9 


ABS*MCD 








0.23±0.1 




2.89 / 2 


6.00E-15 


1.68E+37 


X28=' 


309.8 


ABS*MCD 








0.64±0.1 




24.7 / 17 


3.98E-15 


6.98E+37 


X22 


1411.6 


ABS*PL 


0.6±0.07 




1 Q+0.16 
^^■^-0.14 






86.73 / 85 


4.34E-13 


7.62E+38 



The inner temperature of the multi-eolor disk model. 
^ MEKAL model temperature. 

3 Observed X-ray fluxes and estimated luminosities for the 0.3 — 8 keV energy range, uncorreeted for absorption. 

^ Possible line emission. 

^ Excess flux just over 1 keV. 



Table 4 

HST Images 



HST Instrument Proposal ID Filter Exposure time (s) 

WFPC2 6716 F336W 2x520 

ACS/WFC 10585 F435W 4x3660 (PosA) 

4x3478.91 (PosB) 

F555W 4x2460 

F658N 4x360 

F814W 4x2060 



We also use available F336W (« U) band imaging of 
two pointings within NGC 4449, taken with the WFPC2 
camera (Proposal ID GO: 6716, PI: Stecher). Each posi- 
tion has two exposures. The WFPC2 has four CCDs — 
the Planetary Camera (PC) has a scale of 0."0456 pix~^, 
and the three Wide Field (WF) CCDs have a scale of 
0."0996 pix~^. Note that there are two sets of images 
for each WFPC2 observation available in the HLA: a 
combined WFPC2 image including all four CCDs with 
the PC resampled to the same resolution as the three 
WF CCDs, and an image of only the PC, at its original 
pixel scale. We have used the combined WFPC2 image 
when our cluster candidates were in one of the three WFs 
and WFPC2-PC for objects located in the PC. Details 
on pointings and exposure times for the HST data used 
here are given in Table |4l 

We identified ~200,000 sources in each ACS point- 
ing (in the V band), using the IRAFQ DAOFIND task. 
These include star clusters and bright, individual stars 
in NGC 4449, as well as some foreground stars and back- 
ground galaxies. We perform circular aperture photome- 
try of all detected sources using radii of 1 and 3 pixels and 
background annuli of 8 and 13 pixels, using the PHOT 
task in IRAF. We use zeropo ints on the VEGAMAG 
system taken from Table 9 in iHoltzman et al.l ()1995afl 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



for the WFP C 2/F3 36W filter, and from Table 10 in 
ISirianni et al.l (|2005[ ) for the ACS observations. Cor- 
rections for inefficiency in the charge transfer were de- 
termined for the WFPC2/F336W photometry based on 
the formulae given by Dolphin (2000; we have used the 
most re cent characterization available from A. Dolphin's 
website |http : / /purcell. as . arizona.edu / wfpc2_calib / ) . We 
also applied aperture corrections to extrapolate from our 
fixed aperture radius to the total magnitude in each fil- 
ter, based on the concentration index (C/, the magnitude 
difference measured between a radius of 1 and 3 pixels) 
measured for each source. We obtained a linear fit be- 
tween the measured CI and aperture corrections for «20 
hand-selected, relatively isolated sources, and used this 
fit to estimate the aperture corrections for the rest of the 
sources. 

3.2. Color Magnitude Diagram of Donor Stars 

The colors and luminosities of individual stars give an 
estimate of their masses and hence a constraint on their 
ages. In several cases, the optical images reveal a sin- 
gle point source within the positional uncertainty of an 
XRB, which is likely the donor star. Although a full anal- 
ysis of the spectral energy distributions for these stars is 
beyond the scope of this paper and will be presented in 
a future work, we can use our photometry to broadly 
assess their ages. In Figure IH we compare the mea- 
sured luminosities and colors for these matched optical 
sources with solar metallicity stellar isochrones with ages 
of 10 Myr, 100 Myr, and 1 Gyr f rom the Padova group 
(|Girardi et al.l (f2008h : iMarigoeT aL (2008)). These pro- 
vide a rough guide to the ages/types of the candidate 
donor stars, and suggest that the bright stellar sources 
associated with XIO and X24 are supergiant stars, while 
those associated with X6, X31, and X14 are significantly 
fainter and hence may be somewhat older. 

3.3. Cluster Selection 

We measure t he sizes of all detected objects using the 
Ishape software llLarseni ri999f) . This gives a better mea- 
sure than CI of the sizes of relatively bright, isolated star 
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Figure 3. X-ray color-color diagram of X-ray binaries in NGC 4449. The HI and H2 colors are defined in Section 2. Candidate high 
mass black hole binaries are shown as squares, and the rest of the sample is shown as triangles. The asterisks show the X-ray sources that 
are coincident with clusters. Theoretical tracks are shown as the open symbols connected with solid lines. The green triangles represent a 
power law (PL) with an increasing photon index from 1.0 to 3.0, and the red circles show predictions for diskblack-body models (DBB) with 
increasing temperature from 0.1 to 1.0 keV. The orange squares are a disk black-body model with T=0.9 keV and an increasing hydrogen 
column density (DBB-f-ABS) njj in steps of 0, 1 X 10^°, 5 X 10^°, 1 X 102\ 5 X W'^^ , and 1 X 10^^ cm^^. See text for details. 



clusters. Ishape convolves analytic profiles with the PSF, 
and determines the best fit to each source. We assume a 
KING30 profile, a single mass iKin^ ([l966j model with a 
fixed ratio of tidal to core radius of 30 and fit the data 
within a radius of 5 pixels. For more details on Isha pe 
and estimating the sizes of clusters, see iLarsenl (|1999[ ). 

Most star clusters are slightly broader than the PSF 
at the distance of NGC 4449. The biggest challenge in 
selecting compact star clusters from the entire source list 
is to separate them from chance blends and superposi- 
tions in the crowded star forming regions. We use the 
following critieria to select star cluster candidates: (1) 
my < 22 (My « —6 at the assumed distance of 3.82 Mpc 
for NGC 4449); (2) CI > 1.3; (3) FWHM > 0.2 pix- 
els, i.e. at least 0.2 pixels broader than the PSF, as 
measured by Ishape; and (4) no "neighbors" , i.e. other 



object detected within a 5 pixel radius. A final by-eye in- 
spection was made to throw out some remaining blends 
(at the « 20 — 30% level). The final catalog contains 
129 candidate star clusters, double the nuinber pu blished 
previously bv lGelatt. Hunter &: Gallagheii ()2001|) . Basic 
information for the selected clusters, including their lo- 
cations and measured photometry, is given in Table [5l 
The cluster locations are shown in Figure [2] 

3.4. Cluster Age, Mass and Size Estimates 

Figure [5] compares predictions from solar metallic- 
ity stellar population models of Chariot & Bruzual 
(2007; hereafter CB07, p rivate communication; see also 
iBruzual &: Chariot! (|2003[ )'). with colors measured for our 
final cluster sample. Crosses mark the predicted colors 
for 10^, 10^, 10^, 10^, and 10^° yr clusters, starting from 



XRBs and Star Clusters in NGC 4449 



7 



18 



20 - 



22 - 



24 



26 



1 1 1 1 1 1 

■ 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

^^^^Sl( ftj< 24 1 Myr 


1 1 1 1 1 1 1 1 1 1 1 

■ 


- / 


•X 6 


•X14 


- 






^ 100 Myr 

^ 'XSI / 




- 

: 1 


1 




1 Gyr _ 






Jl 


< < 1 







1 2 

(l^FSSSW " n^F814w) 


3 



Figure 4. Color magnitude diagram of optical point sources which 
are unique counterparts to X-ray sources. The solid lines are 
Padova isochrones (for Z = 0.02) of 10 Myr, 100 Myr and 1 Gyr 
(from left to right). 



Table 5 

Star cluster catalogue 



Number 


RA 


DEC 


U-B 


V-I 


Age^ 










1 


3g(T/yr) 


1 


187.05396 


44.083114 


-2.106 


-0.295 


6.40 


2 


187.05370 


44.083248 


-2.195 


-0.395 


6.02 


3 


187.06067 


44.083354 


-0.303 


1.080 


9.30 


Note. — 


This table 


is published 


in its entirety in the 


electronic 



edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

^ Typical uncertainties arc approximately ±0.3 in log r. 



the upper left. The arrow shows the d irection of redden - 
ing for a Galactic- type extinction law ([Fitzpatrickll 1 9991 ) . 

The model predictions match our cluster photometry 
relatively well. The majority of the clusters have blue 
integrated colors, indicating that they are fairly young, 
with ages r <^ few x 10^ yr. At the top of the dia- 
gram there are a few young clusters in crowded regions 
that fall above the models, likely due to some contam- 
ination from neighbors in the lower resolution U band 
images ( we see a simila r effect for the colors of clusters 
in M51; IChandar erahl (|2010aD ). The few sources with 
U-B colors significantly redder than the model predic- 
tions are somewhat faint in the U band, leading to large 
uncertainties in this color. The bottom panel of Fig- 
ure [5] shows the B — V vs. V — I two-color diagram, 
and includes some objects for which no U band photom- 
etry was measured, either because they are outside the 
field of view or because they are too faint in this band. 
The sources which fall significantly below (redward of) 
the model predictions in B — V have strong nebular line 
emission (which is observed in the narrow-band F658N 
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Figure 5. Optical two-color diagrams of all star clusters in our 
sample. The solid line shows predictions from the single stellar 
population models of CB07 and the crosses show the following 
predicted ages: log(T/yr) = 6, 7, 8, 9 and 10, starting from the 
upper-left. The arrow shows the direction of reddening. 



(Ha) filter), moving them off the model track. While 
two-color diagrams like those shown in Figure [5] are use- 
ful for visualizing the evolution of clusters, we actually 
use a minimization technique to estimate the ages of 
the clusters. This means that points outside the models 
(e.g., the "high" points in Figure [5]) can be reasonably 
well fit. 

We estimate the age r and extinction Ay for each clus- 
ter as we have done in previous works (see e.g.. lFall et al.l 
(^005) for details), by performing a least fit comparing 
observed magnitudes with the predictions from the CB07 
stellar population models with solar metallicity {Z — 
0.02), which appear to match the measured colors of clus- 
ters in NGC 4449 reasonably well. The best fit combina- 
tion of T and Av for each cluster returns the minimum 

X^: X^{r,Av) = Y^x^y^i^^f"^ '^T^^f ^ where the sum 
runs over all available broad-band (UBVI) filters and the 
F658N narrow-band filter, but requires a minimum of 
three measurements (including the V band) to estimate 
age and extinction. The weights W are related to the 
photometric uncertainty a\ as W\ = [cr| -I- (0.05)^]"^. 
The F658N filter includes both stellar continuum and 
nebular line emission, and is dominated by line emission 
from ionized gas for the youngest (t ^ several x 10® Myr) 
objects, and by continuum emission from stars for clus- 
ters with ages of r ^ 10^ yr. This enables us to use the 
narrow-band filter as a fifth data point in many cases, 
regardless of the age of the cluster. 

We tested our method by estimating ages and masses 
from our analysis for different assumptions regarding 
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extinction and by comparing our photometry with the 
predictions from a model with lower metallicity {Z — 
0.008). In general, we find that the age estimates are 
quite robust for very young clusters (t ^ 6 x 10^ yr) 
that have nebular emission, and for clusters with ages 
between 3 x 10^ yr<;, r ^ 10^ yr. Clusters older than 

« 10^ yr suffer from the well-known age-metallicity de- 
generacy, and we cannot tell if they are approximately a 
Gyr or older than « 10 Gyr from the data presented here. 
Eleven clusters in our catalog have integrated colors simi- 
lar to those of metal-poor globular clusters in our Galaxy 
(and in many other galaxies) with 0.8 <V — I < 1.3 and 
0.6 < B — V < 0.9. The dating of clusters with no neb- 
ular emission and with ages between w 6 x 10^ yr and 
« 3 X 10^ yr is degenerate in some cases, with two nearly 
equally good combinations of age and extinction. We es- 
timate that typical uncertainties are w 0.3 in log t for 
clusters with r ^ 10^ yr, typical for this method (e.g., 
iFall et all ^00^ . 

The mass of each cluster is estimated from the total V 
band luminosity, corrected for extinction, and the age- 
dependent mass-to-light ratios [M / Ly) predicted by the 
CB07 models. 

4. FORMATION AND DISRUPTION OF THE 
CLUSTERS 

The mass and age distributions of a population of star 
clusters provides important clues to the formation and 
disruption of the clusters. In Figure [6] we show the mass 
and age estimates for star clusters in NGC 4449. The 
solid line represents My = — 7, the approximate com- 
pleteness limit of our sample, and shows that our sample 
does not contain clusters over the same mass range at all 
ages, because clusters fade over time. The mass-age dia- 
gram shows some small scale features, such as a sparsely 
populated region in the range 7.0 <^ log(r/yr) <^ 7.6 yr. 
This particular feature occurs where the predicted colors 
loop back on themselves, covering a small region in color 
space over a relatively long time, and effectively result- 
ing in a gap. This artificial, empty stripe and similar 
features do not affect our conclusions. 

Qualitative trends in the distribution of cluster ages 
and masses are apparent from the mass-age diagram. We 
see that NGC 4449 has formed relatively massive clus- 
ters {M J> 10^ Mq) more or less continuously over the 
last 10^ yr. A number of clusters (approximately half of 
our sample), formed very recently, in the last t <^ 10^ yr, 
and a number of clusters clearly formed w 1 — 3 x 10^ yr 
ago. As mentioned in Section 3.3, our sample includes 
« 11 candidate ancient globular clusters, implying that 
NGC 4449 began forming stars approximately a Hub- 
ble tim e ago. This is broadly consistent with the re- 
sults of lAnnibali et al.l ()2008f ). who found that stars in 
NGC 4449 have formed over at least the last « 10^ yr, 
with tentative evidence for earlier star formation as well, 
based on the colors and luminosities of individual stars. 
Qualitatively, the overall distribution of cluster ages and 
masses in NGC 4449 appears to be similar to that in 
other galaxies with very different environments and star 
formatio n histories , such as the merging Antennae galax- 
ies (e.g.. iFall et al.' (2005*), and the more quiescent Mag- 
ellanic Clouds ()Chandar et al... 2010a) . although the gap 
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Figure 6. Mass-age diagram of all clusters in our sample. The 
solid line shows the approximate magnitude limit of My = —7.0. 



between 10*" — 10® yr is somewhat more prominent. We 
present a more quantitative analysis below. 

Figure [7^ shows the age distribution for clusters in 
NGC 4449 in two different intervals of mass. Note that 
the data points for each distribution are restricted to 
a mass-age range that is not affected by incomplete- 
ness (i.e. above the solid line), so the observed shapes 
are caused by the formation and disruption of the clus- 
ters and not by their fading out of our sample. Al- 
though there are few data points, these distributions ap- 
pear to decline steeply, and can be approximated by a 
simple power law, dN/dr oc t"' , with best fits of 7 = 
—0.91 ±0.13 for clusters more massive than 2 x lO'* Mq, 
and 7 — —1.21 ± 0.21 for clusters with masses between 
1 — 2x10"* Mq . These values of 7 are the same within the 
uncertainties. We find that 7 can change by up to ~ 0.3 
if different bin sizes and centers are used, but the over- 
all decline starting at very young ages does not change. 
Based on these experiments, we take 7 = —1.0 ± 0.3. 

Figure [TId shows the mass function for clusters in 
NGC 4449 in two different intervals of age. These 
mass functions can be approximated by a power law, 
dN/dM oc M^, with best fit values of ^ = -2.25 ± 0.15 
for T < 10^ yr, and /3 = -2.07±0.19 for r = l-2x 10® yr, 
which are the same within the uncertainties. We take /3 
to be — 2.16±0.3, the mean and range of the values above. 

Because the mass and age distributions appear to be 
at least approximately independent of one another, the 
joint distribution of cluster ages and masses g{M, r) can 
be approximated as g{M,T) oc M^t'^, with P sa —2 and 
7 « — 1, for r ^ few X 10® yr and M ^ 10^ Mq, similar to 
the form observed for cluster s in the Magellanic C louds 
(jChandar et al.l2010aD. M83 (iCha ndar et al.l20T0bl). an d 
i n the Antennae (iWhitinore et al . 2007; F aff et al.ll20Tol) . 

IC handar et al.l ( 2010cD compiled results for the age dis- 
tribution of star clusters in over a dozen different galax- 
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Figure 7. The age distribution of clusters in NGC 4449 in the indicated intervals of mass is shown on the left, and the mass function in 
the indicated intervals of age is shown on the right. The lines show the best fits to the distributions. See text for more details. 



ies, including dwarf irregular, spiral, and merging galax- 
ies. All of these different galaxies appear to have similar 
shapes to their (mass-limited) age distributions, where 
they decline from the present to the past (over the past 
« fewx 10^ — 10^ yr). This shape for the age distribution 
has previously been interpreted as due primarily to the 
gradual, early disruption rather than to the formation of 
the clusters, since it is far more likely that the clusters 
in all of these galaxies have similar disruption histories 
than it is that they have similar formation histories. 

After their formation in the dense cores of giant molec- 
ular clouds, different physical processe s cause clusters to 
lose mass and to eventual ly disrupt. iFall et all ()2009l ) 
and lChandar et ahl (|2010al ) have suggested the following 
approximate sequence and timescale for these processes: 

(1) removal of ISM by stellar feedback, r < 10^ yr; 

(2) continued stellar mass loss, 10^ yr ^ t < 10* yr; 
and (3) tidal disturbances by passing molecular clouds, 
r > 10* yr. It is likely that the cluster mass and age, 
which can be approximated by power-law distributions, 
result from a complex situation that involves several of 
th ese disrupt ion process es. We refer the interest ed reader 
to IFall et al.l (pOOal and lChandar et all (|2010aD for more 
details. On longer timescales, mass-loss is driven by the 
escape of stars due to internal two-body relaxation, or 
evaporation. 

5. SPATIAL CORRELATION BETWEEN X-RAY 
AND OPTICAL SOURCES 

5.1. Correlation Between the Positions of XRBs and 
Star Clusters 

iKaaret et al.l (|2004f ) previously suggested that X-ray 
sources in three different starburst galaxies (M82, NGC 
1569 and NGC 5253) may have formed in young star 
clusters, because they are preferentially located near star 
clusters, albeit with a significant displacement (« 200 pc 
on average). This conclusion was based on a compari- 
son between the cumulative distribution of displacements 
measured between actual X-ray sources and their clos- 



est clusters, and simulated displacements from a popu- 
lation of X-ray sources distributed randomly relative to 
the clusters. In order to explore whether or not X-ray 
binaries form in compact star clusters, here we compare 
the positions (of the entire sample) of 23 candidate X-ray 
binaries (presented in Section 2) with those of the star 
clusters (Section 3), to look for a statistically significant 
spatial correlation between the two. 

We repeat the spat ial correlation analysis performed 
bv IKaaret et~all (|2004l ) for our sources in NGC 4449. We 
identify the nearest star cluster to each X-ray binary, and 
measure the (projected) distance between the two. For 
XRBs that are not coincident with clusters, this gives a 
lower bound on the displacement of the source from its 
parent cluster, since the closest cluster is not necessar- 
ily the parent. The results for NGC 4449 are shown in 
Figure [8] as a cumulative distribution of the separations 
between the binaries and the closest cluster. The figure 
shows that there are 7 XRBs that are within 100 pc of a 
star cluster, and nearly all are within w 400 pc of a star 
cluster. Interestingly, three HMXBs in NGC 4449 are 
spatially coincident with very young star clusters (within 
the la positional uncertainty of « 0.5" for Chandra), 
suggesting that these XRBs formed within the clusters. 

To determine whether the clustering of XRBs near star 
clusters is statistically significant, we compare the obser- 
vations with simulated X-ray source populations. We 
consider two different randomly distributed cases, i.e. 
where the X-ray sources are not associated with the clus- 
ters. First, we generate rand om sets of (23) unif ormly 
distributed sources, as done in IKaaret et al.l (|2004f ). and 
find spatial displacements from the star clusters using 
the procedure described above. The simulation was run 
1000 times, with the cumulative distribution determined 
in exactly the same way as for the actual observations, 
and the median result shown as the blue line in Figure |8l 
The uniformly distributed random sources show a clear 
difference from the observ ed spatial disp l aceme nts, simi- 
lar to the results found bv IKaaret et al.l (|2004D for three 
different starburst galaxies (not including NGC 4449). 
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Figure 8. Cumulative distribution of the displacement (in par- 
sees) between X-ray binaries and star clusters. The black (solid) 
line shows the result for sources in NGC 4449. The other lines show 
results from our Monte Carlo simulations. The red solid line is the 
median of 1000 random distributions where the the locations of the 
X-ray sources are drawn from a distribution that follows the light 
of the galaxy. The dashed lines show the 1 a deviation from this 
median curve. The blue solid line represents the median of 1000 
random distributions where the X-ray source locations are drawn 
randomly from across the face of the galaxy. 

Next, we select random sets of (23) sources distributed 
in a more realistic fashion, one where the XRBs fall off 
radially like the stellar light in NGC 4449 itself. This 
was accomplished by fitting the radial luminosity profile 
of NGC 4449 in a 2MASS K band image, by a broken 
power-law, 4>{L) oc L", with a = —0.79 over the radial 
range 150 — 800 pc, and a — —2.28 over the radial range 
800 — 3000 pc. The K band should give a reasonably 
good measure of the distribution of older stars that form 
the stellar backbone of NGC 4449, rather than the dis- 
tribution of young star clusters. The median result from 
1000 runs of the random simulations described above are 
shown as the solid red line in Figure |S1 In this case, the 
results of the randomly drawn spatial displacements are 
more similar to the observations, except at the small dis- 
placement end, with the observations now overlapping 
partly with the la uncertainty line from the simulation 
at separations > 100 pc. We find that a radial profile of 
NGC 4449 determined from far-ultraviolet images taken 
with GALEX gives very similar results to those based on 
the K-band profile. Our results indicate that the spe- 
cific recipe used to populate synthetic X-ray sources in 
the image can affect the interpretation of the results, 
i.e. whether or not the X-ray sources are associated with 
star clusters. In general, a uniform distribution is much 
more likely to place synthetic X-ray sources in the outer 
portions of the galaxy when compared with the second 
method; because there are few X-ray sources located in 
the outskirts of NGC 4449, this procedure causes the 
simulation and observations to look substantially differ- 
ent. 



We conclude that, in the case where the sources are 
distributed randomly but in a fashion similar to the un- 
derlying galaxy light, the random distribution does not 
differ significantly from the observed one, except at the 
smallest separations, where we find several XRBs closer 
to star clusters than predicted by the random distribu- 
tions. These sources are discussed in more detail in Sec- 
tion 6. 

5.2. Optical Sources That are Coincident With 
HMXBs 

Three of the candidate HMXBs (X12, X21 and X46) 
have young (r <^ 10^ yr) star clusters located within their 
« 0.5" positional uncertainties. There is also an XRB co- 
incident with an old globular cluster (X29); this is almost 
certainly a low-mass X-ray binary (LMXB). Here, we 
check the probability that this spatial coincidence is due 
to chance superposition rather than to a physical associa- 
tion between the HMXBs and the clusters. We distribute 
23 XRBs (the number in our sample) randomly through- 
out NGC 4449 in 100,000 runs, and find one chance su- 
perposition approximately every fifth run when all 129 
clusters were included. This is «25 times less than the 
four real observed coincidences between star clusters and 
XRBs. The probability of a chance superposition is al- 
most identical if we only consider the three coincident 
HMXBs and star clusters younger than 10 Myr, i.e. the 
most probable host clusters. Here, the random sources 
follow the luminosity profile of NGC 4449 (see Section 
4.1). If the random sources are uniformly distributed 
instead, we find coincidences only ^3% of the time, ap- 
proximately 1 out of every 30 runs. Our results indicate 
that the positional coincidence between the XRBs and 
clusters in these three cases has a high statistical sig- 
nificance, and it is highly unlikely the result of chance 
superposition. 

In addition to the XRBs that are coincident with star 
clusters, several others are coincident with an optical 
point source: X6, XIO, X14, X20, and X24, as presented 
in Section 3.2. These point sources are well within the 
body of NGC 4449. They have the colors and magnitudes 
expected of supergiant stars, and are almost certainly 
the donor star in the X-ray emitting binary. Other can- 
didate HMXBs are located in crowded regions, making it 
more difficult to identify a unique optical counterpart to 
the XRB, because there are several sources in the HST 
image within the positional uncertainties of the X-ray 
source. Some of the X-ray sources (e.g., X18, X13, X20, 
XI 1, and X19) do not have any obvious optical counter- 
part. 

6. DISCUSSION 

6.1. Properties of Star Clusters Closest to HMXBs 

The ages of star clusters in the vicinity of the X-ray 
binaries can help to constrain the ages and the nature of 
the HMXBs. Even when there is no unambiguous stellar 
or cluster counterpart to an X-ray source, it is plausible 
that the XRB is the same age as the stars and clusters 
around it. Here, we present the physical properties of 
the star cluster that is closest to each HMXB. In Table 
[HI for each HMXB we compile the distance to the closest 
star cluster in our catalog, and the age and mass of this 
cluster. We also include the X-ray luminosity of each 
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source determined in Section 2.1, and a column that lists 
whether the X-ray colors are better described by a disk 
black body or a power law model, or are in the absorbed 
portion of the diagram. 

The colors of HMXBs can vary, because these binaries 
go through different emission states, or because of their 
geometric orientation, since inclination effects through 
a disk can result in colors in the absorbed portion of 
the X-ray two-color diagram. This means that it is 
not possible to definitively identify the type of HMXB 
(e.g., black hole vs. neutron star) or to determine its 
age, just from the X-ray properties. However, different 
types of HMXBs do show some general trends in the 
color-color diagram. For example, black hole binaries 
often approximately follow the disk black-body models 
(jRemillard fc McClintockll2006l ). At high accretion rates 
the X-ray emission from black hole binaries is dominated 
by thermal emission from a disk, leading to both high 
luminosities and soft X-ray colors. On the other hand, 
accreting low-mass neutron star binaries often have ab- 
sorbed power-law spectra. 

6.2. A Population of Very Young, Massive Black Hole 
Binaries in NGC 444^ 

We can select very young binaries, which are likely to 
be high-mass black hole binaries (BHB), in two different 
ways: (1) from their location in the X-ray color-color 
diagram, and (2) from their proximity to very young 
star clusters. We first select candidate BHBs from Table 
[6] based on their X-ray colors. There are seven X-ray 
sources that are best described by disk black-body mod- 
els: X12, X21, X23, X24, X28, X42 and X46. All of these 
X-ray sources are either coincident with (X12, X21, and 
X46), or fairly close to (within 200 pc) a very young 
r ^ 8 Myr star cluster^ The coincidence with clusters 
is particularly important, because it establishes a direct 
connection between the BHBs and very young clusters 
(recall that in Section 5.1 we found that chance superpo- 
sitions are highly unlikely). These seven sources have a 
median X-ray luminosity of 2.25 x 10^^ erg s~^. The fact 
that they have the X-ray colors expected for BHBs and 
are also close to very young star clusters strongly sup- 
ports that these sources are in fact, very young BHBs 
(discussed in more detail below). 

Next, we select candidate black hole binaries based 
solely on their proximity (within 200 pc) of a very young 
(r ^ 8 Myr) star cluster. These criteria return all seven 
XRBs described by the disk black body models, plus four 
additional sources: X6, X22, X43, and X44. All four 
of these are found in the "absorbed" portion of the X- 
ray two-color diagram. However, they are located near 
star forming regions which contain very young t <^ 8 
Myr star clusters, and that are further out in NGC 4449 
(with galactocentric distances Rg ^ 1.2 kpc) on aver- 
age than the seven sources discussed previously. The 
stellar density of the galaxy has dropped significantly at 
these galactocentric distances, which strongly suggests 
that X6, X22, X43, and X44 are associated with the 
nearby, recent star formation. In fact, random simula- 
tions return only a single XRB with Rg ^ 1.2 kpc and 

* Source X21 also happens to have an older, t 200 Myr cluster 
nearby. 



within 200 pc of a young star cluster in 10,000 runs. 

The physical association between BHBs and very 
young star clusters is also apparent as the small sepa- 
rations in the cumulative distribution of spatial displace- 
ments between XRBs and star clusters (Figure ^ . This 
regime deviates strongly from random distributions, as 
we would expect if BHBs are physically associated with 
very young star clusters. We conclude that (at least 
some) BHBs form in, and not just near, compact star 
clusters, based on the fact that three candidate BHBs 
in our sample are spatially coincident with a very young 
star clusters and that a fourth is very close to a young 
cluster (within a projected distance of 13 pc). 

The ages of r <^ 6 — 8 Myr estimated for the coincident 
and proximate star clusters to the 11 XRBs discussed 
above strongly suggests that these are black hole (rather 
than neutron star) binaries. At these very young ages, 
only stars initially more massive than M ^ 25 — 30 Mq 
will have had time to become supernovae (based on the 
Padova models for solar metallicity) . While there is still 
some uncertainty about the exact range of stellar masses 
that end their lives as black holes rather than as neutron 
stars, most models predict that the transition between 
these two types of compact remnants occurs for star with 
initial masses somewhere in the range of f» 18 — 25 Mq 
(e.g.. Fryer (1999)), below the initial stellar masses that 
have completed their main sequence lifetime in ss 6 — 
8 Myr old star clusters. However, metallicity effects may 
complicate the relatio nship between initi al stellar mass 
and remnant type (e.g.. lHeger et al.l ()2003D ). Given these 
uncertainties we will refer to these as 'candidate' BHBs. 

One of the XRBs that is best described by a disk black- 
body model, X24, is not coincident with a star cluster, 
but is coincident with a bright point source. The lack 
of any other bright stars within the la astrometric un- 
certainty suggests that this is the high mass donor star 
in the XRB. The luminosity and color of this source are 
consistent with isochrones that are « 8 — 10 Myr, as 
shown in Figure IH similar to the age of the nearest star 
cluster. 

To summarize, the main results of this section are that 
we find strong evidence for a population of BHBs in 
NGC 4449, and that these massive binaries likely formed 
recently in compact star clusters. Three of the can- 
didate BHBs appear to reside within their parent star 
clusters, while the others do not. These very young 
(r <^ 6 — 8 Myr) sources comprise a significant fraction 
of all X-ray emitting binaries brighter than w few x 10^^ 
erg s~^ in NGC 4449. BHBs therefore make up approx- 
imately 48% or 11 out of 23 XRBs. 

6.3. Processes Responsible for the Spatial 
Displacement Between BHBs and Star Clusters 

Black holes are the compact remnants of massive O 
stars. O stars typically form in (massive) clusters rather 
than in the field, although they can be dynamically 
ejected from clusters into the field. For e xample, stud- 
ies o f field O stars in the Milky Way (|de Wit et al.l 
120051 ) show that the majority of field O stars appear 
to be runaways from nearby star clusters. Moreover, 
many and possibly most O stars are bo rn in binaries 
(|Garcia fc MermiliiodI [2OOII: iLarsonl [200l . Section 6.2 
provided strong evidence that BHBs probably form in 
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Table 6 

Properties of Star Cluster Closest to XRBs in NGC 4449 



ID 


X-ray 
Model 


ID of Closest 
Cluster 


Cluster Age-"- 
log(T/yr) 


Distance 
(pc) 


Notes 


X8 


Absorption 


96 


9.34 


266 




X14 


Absorption 


117 


8.06 


339 




X18 




117 


8.06 


285 




X13 


Absorption 


112 


8.06 


232 




X6 


Absorption 


93 


6.64 


122 




X29 




89 


8.91 


7 


LMXB 


XIO 


Power-law 


99 


8.31 


286 




X12 


Disk blackbody 


104 


6.96 


10 


Coincident 


X20 


Absorption 


33 


8.86 


222 




Xll 


Absorption 


102 


8.36 


89 




X23 


Disk blackbody 


66 


6.58 


181 




Vic; 
AID 


Power law 


lo 


8.31 


9 


Coincident 


X24 


Disk blackbody 


69 


6.84 


68 




X21 


Disk blackbody 


53 


6.72 


13 


Coincident 


X28 


Disk blackbody 


81 


6.72 


136 




X31 




87 


6.70 


630 




X22 


Absorption 


70 


6.38 


133 




X19 


Super soft source 


51 


8.06 


383 




X41 


Super soft source 


50 


6.90 


225 




X42 


Disk blackbody 


6 


6.40 


113 




X43 


Absorption 


62 


6.40 


57 
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Uncertainties in the cluster age estimates are typically ±0.3 in log r. 

star clusters. 

In the previous section, we found that while three of 
the candidate BHBs in NGC 4449 are coincident with a 
star cluster, the majority of them are not. What physical 
process(es) is responsible for this displacement between 
BHBs and their parent clusters? We consider three dif- 
ferent mechanisms that could lead to the apparent dis- 
placement between BHBs and their parent star clusters: 

• The parent cluster has dissolved and is therefore no 
longer visible, leading to an apparent displacement 
between the BHB and neighboring clusters. 

• The BHB was ejected from its parent cluster during 
dynamic al interactions with stars in the dense clus - 
ter core (|Poveda et al.l[T967t iGies fc Bolton|[T986h . 

• The BHB was ejected from its pare nt cluster due 
to an asymm etric supernova kick (jZwickvl 119571 : 
IBlaauwIlTgel . 

In Section 4, we suggested that the shape of the clus- 
ter mass and age distributions in NGC 4449 are primarily 
signatures of the disruption rather than the formation of 
the clusters. Regardless of whether individual clusters 
dissolve partially or fully, our results in NGC 4449 (and 
those in other galaxies), suggest that star clusters lose a 
significant amount of mass very quickly, on timescales of 
only « 10 Myr0 UV spectra of starburst galaxies reveal 
that the dispersed field populations are dominated by 
B stars, whereas UV- bright star clusters are often dom- 
inated by O stars, consistent with a scenario of rapid 

^ This early mass loss is not driven by relaxation of t he clu s- 
ter due to two-body interactions (e.g., IFall & Zhanel I I2001I) ). 
Relaxation-driven evaporation operates on significantly longer 
timescales than the ages of the BIIBs, and ejects mostly low mass 
stars. 



cluste r dispersal and / or mass-loss (e.g., iTremonti et al.l 
(pOOl[ ): [ChMldar et al.l (2003)). Despite growing evidence 
that even massive star clusters, i.e., those most likely to 
host O stars, may disperse rapidly, we believe that this is 
unlikely to be the mechanism responsible for the observed 
spatial displacement between very young BHBs and star 
clusters. N-body simulations show that an unbound clus- 
ter retains the appearance of a bound cluste r for 10 — 20 
crossing times ([Baumgardt &: Kroupal [20071) . on the or- 
der of 10-20 Myr for typical clusters in our NGC 4449 
catalog. Since no cluster is observed at the locations of 
the majority of the BHBs, the early dispersal of clusters 
can be ruled out as the origin of the spatial displacement 
between very young, high mass BHBs and star clusters. 

We conclude that many BHBs in NGC 4449 have been 
ejected from their parent clusters, either via dynamical 
kicks due to interactions with other stars in cluster cores, 
or due to an asymmetry in the supernova explosion of the 
compact object, which nonetheless does not unbind the 
binary. We cannot differentiate between these two mech- 
anisms using the current observations. We can however, 
estimate lower and upper limits to the ejection velocities. 
We estimate a lower kick velocity by assuming that each 
BHB was ejected from its closest cluster very soon after 
it formed, and divide the distance to this cluster by its 
age, neglecting for the moment uncertainties in the clus- 
ter age estimates. These lead to lower limits of 45, 10, 
27, and 23 km/s for X23, X24, X28, and X42, respec- 
tively. The other BHB binaries give limits of 28, 55, 23, 
and 84 km/s for X6, X22, X43, and X44. The cluster 
ages and hence the velocities are uncertain by w a factor 
of two, resulting in lower limits to the velocities between 
« 5 — 160 km/s. Of course dynamical ejection need not 
have occurred right after the cluster formed, but could 
have occurred more recently. We estimate upper limits 
to the kick velocities by assuming that each BHB was 
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ejected within the last 1 Myr, i.e. we divide the distance 
by 1 Myr. This procedure gives a range of velocities (for 
non-coincident sources) between « 30 — 180 km/s. 

6.4. The Nature of Older X-ray Binaries in NGC 4449 

In the previous subsections we found that w 11 of the 
23 candidate XRBs in NGC 4449 are likely very young, 
massive BHBs. The goal of this section is to better un- 
derstand the ages and types of the remaining X-ray bi- 
naries. 

We first constrain the number of low mass X-ray bina- 
ries (LMXBs) in NGC 4449, systems consisting of a black 
hole or neutron star accreting from a low-mass compan- 
ion. Source X29 is almost certainly a LMXB, because it is 
spatially coincident with a cluster that has integrated col- 
ors similar to those of ancient Galactic globular clusters. 
Its X-ray colors and luminosity however, are unremark- 
able when compared with the rest of the XRB sample, 
and therefore cannot be used to select ancient LMXBs in 
general. In their study of five early-type galaxies (ellip- 
ticals and len ticulars) , whi c h hav e formed LMXBs but 
not HMXBs, IKundu et al.) (l2007l) find « 25 - 60% of 
all LMXBs are coincident with globular clusters and the 
rest are found in the field. Assuming a similar fraction of 
cluster-to-ficld LMXBs in NGC 4449 implies that there 
are « 2 — 4 field and hence « 3 — 5 total LMXBs in 
NGC 4449, « 13 - 22% of our total sample. 

Our analysis implies that there are 7 — 9 remaining X- 
ray binaries in NGC 4449 that are neither very young, 
high mass BHBs nor very old LMXBs. These are likely 
older than w 10 Myr and younger than several billion 
years (age of LMXBs). Interestingly, from Table [6] we 
find that the closest clusters to X14, X18, X13, XIO, 
X20, XI 1, and X19 are all intermediate age, with r w 
100 — 400 Myr. The separations between these X-ray 
binaries and their closest clusters is larger than those 
found between BHBs and their closest clusters, with a 
median (mean) separation of « 285 pc (250 pc). 

These "intermediate" age HMXBs are likely neutron 
star X-ray binaries with either a Be or supergiant star 
for a companion. Accretion from the former occurs 
via a disk, while the latter are primarily wind fed. 
There are also differences in the optical luminosities 
of Be versus supergiant donor stars: Be-HMXBs in 
the Sm all Magellanic Cloud have My ranging from —2 
to -5 (jMcBride et all 120081 ). which at the distance of 
NGC 4449 corresponds to my ~ 26 — 23, while super- 
giant HMXBs in the Galax y are significantly brighter 
(jChevalier fc Ilovaiskvl |1998[ ), with My brighter than 
« —6.5 (my ^ 21.5). XIO and X14 are therefore can- 
didate supergiant-HMXBs, because they have bright co- 
incident point sources (see Section 3.2). X18, X13, X20, 
XI 1, and X19 however, do not show any bright coincident 
point sources in the optical images down to my ~ 25, 
and are therefore likely Be-HMXBs. 

Be- X-ray binaries domin ate the population of HMXBs 
in the Magellanic Clouds (iMcBride et~all I2008D . They 
turn-on approximately 20-50 Myr after a star formation 
event, as the first neutron stars are formed. The low mass 
of the neutron star and lower accretion rates mean that 
Be- X-ray binaries typically have lower X-ray luminosi- 
ties than massive BHBs. However, most Be star X-ray 
binaries go into outburst when the neutron star passes 
through the circumstellar disk of the Be star and the 
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Figure 9. Same as Figure [S] but excluding the eleven massive 
BHBs and the known LMXB. 



accretion rate increases. We note that the candidate Be- 
HMXBs in NGC 4449 have X-ray luminosities consider- 
ably higher than quiescent Be X-ray binaries. If they 
are indeed Be X-ray binaries, then they are in outburst, 
and we are detecting only the highest luminosity sources. 
This would suggest there is a much larger population 
of Be X-ray binaries in NGC 4449 below our detection 
threshold. 

In Section 6.2 we used the spatial displacement dia- 
gram to support our conclusion that very young, massive 
BHBs form in compact star clusters. Here, we note that 
there is an inherent limitation in our ability to interpret 
this diagram for older binaries. While the association 
between very young clusters and XRBs shows up clearly, 
this is not true for intermediate age sources. Figure 9 
shows the cumulative distribution of displacements be- 
tween T ^ 10 Myr X-ray binaries and star clusters, where 
we have removed the BHBs and the known LMXB. After 
w 50 — 100 Myr, XRBs that have been ejected at even low 
velocities of a few km/s are no longer obviously spatially 
associated with their parent clusters, since at separations 
larger than « 200 pc the distribution is consistent with 
the 1(7 uncertainties from a random distribution. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we presented the discovery of 23 can- 
didate X-ray binaries in the nearby starburst galaxy 
NGC 4449, from C/iandra/ ACIS-S observations. We 
measured count rates, luminosities, and colors for these 
sources. 

We also presented a new catalog of 129 compact star 
clusters brighter than My « -7 in NGC 4449 from 
multi-band, high resolution, optical imaging taken with 
the ACS/WFC and WFPC2 cameras on-board the HST. 
This doubled the number of compact star clusters known 
in this galaxy. Mass and age estimates for these clus- 
ters show that NGC 4449 has formed relatively massive 
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(M ^ IO^Mq) clusters more or less continuously over 
the last w 10^ yr. The joint distribution of cluster ages 
and masses appears to be similar to those found in other 
nearby galaxies such as the Magellanic Clouds, M83, and 
the merging Antennae galaxies, albeit with somewhat 
larger uncertainties due to the smaller number of clus- 
ters, and can be approximated as g{M,T) oc M^t^ , with 
^ = -2.16 ± 0.30 and 7 = -1.0 ± 0.3. 

Our main conclusion is that we have found clear evi- 
dence for a population of very young, high mass, black 
hole X-ray binaries in NGC 4449. We find 11 candi- 
date high mass BHBs, nearly half of the sample of X-ray 
emitting binaries. Three of the BHB candidates arc coin- 
cident, within the astrometric uncertainties of Chandra, 
with a very young, r ^ 8 Myr star cluster, and a fourth 
is nearly coincident with a very young cluster. The oth- 
ers are all within 200 pc of very young clusters. Based 
on these results, we suggest that these massive BHBs 
form in star clusters, where most massive O stars (the 
progenitors of BHBs) are born. Many are subsequently 
ejected from their parent clusters either due to dynami- 
cal interactions within dense clusters or as the result of 
an asymmetric supernova explosion. The observed dis- 
placement between BHBs and very young star clusters 
is not caused by the dissolution of the parent clusters. 
The small separation between massive BHB candidates 
and their closest cluster clearly deviates from randomly 
populated distributions, further supporting our conclu- 
sion that these X-ray binaries have a direct relationship 
with young star clusters. 

We found one X-ray binary in NGC 4449 that is coin- 
cident with an old star cluster, and hence is almost cer- 
tainly a LMXB. Based on the fraction of LMXBs found 
in the field vs. in globular clusters in early-type galax- 
ies, we estimate that there are w 2 — 4 field LMXBs and 
hence 3-5 total LMXBs in NGC 4449, although they 
cannot be identified from their X-ray properties alone. 

We suggest that the remaining XRBs are intermediate 
age supergiant and Be-HMXBs. Although none of these 
remaining binaries arc coincident with a star cluster, the 
closest cluster to these sources has an estimated age be- 
tween 100 — 400 Myr. In two cases, we identify coincident 
point sources in the HST images that have luminosities 
consistent with those expected of supergiant donor stars. 
The other candidate intermediate-age X-ray binaries do 
not have an obvious optical counterpart, consistent with 
the expected luminosities for Be star donors. 

Based on these results, we conclude that high mass 
X-ray binaries, particularly massive, black-hole binaries, 
dominate the X-ray luminosity from NGC 4449. We sug- 
gest that the ages and locations of star clusters provide 
important insight and constraints on the different types 
of X-ray binaries in nearby star-forming galaxies. 

We thank the anonymous referee, whose careful read- 
ing and helpful suggestions significantly improved our 
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